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Abstract 

We investigate the transverse spatial profile of down-converted light 
produced by noncollinear, degenerate, Type-I spontaneous parametric 
down-conversion in two types of nonlinear crystals. We hnd that the 
pattern produced by one crystal, beta barium borate (BBO), produces 
a circular down-conversion pattern while the other crystal, bismuth tri¬ 
borate (BiBO) produces an elliptical pattern. We show this difference is 
due to the angle-independent refractive index experienced by the daugh¬ 
ter photons in BBO, while they experience an angle-dependent refractive 
index in BiBO. We image the transverse spatial prohle of the generated 
light to determine the eccentricity produced by each crystal and develop a 
model to explain our observation. Among other things, this model predicts 
that there is a wavelength for which the eccentricity from BiBO is nearly 
zero. Finally, we discuss how the elliptical ring pattern produced in BiBO 
potentially affects polarization entanglement for experimental setups that 
collect biphotons around the entire down-conversion ring. We show that 
the quality of polarization entanglement as measured by the overlap inte¬ 
gral of the spectrum of the two rings, can remain high (> 99.4%) around 
the entire ring at the expense of decreased biphoton rate. 


1 Introduction 

The nonlinear optical process of spontaneous parametric down-conversion (SPDC) 
is a popular method for producing correlated photon pairs, also referred to as 
biphotons or daughter photons. A major advantage of this process is the ability 
to quantum mechanically entangle the daughter photons in various degrees- 
of-freedom including polarization, time, frequency, space, and momentum, for 
example. The entanglement quality can be very high, making it an attrac¬ 
tive source for applications in quantum communication, such as quantum key 
distribution (QKD) [HIS [3], as well as experiments in fundamental quantum 
information science 13 HI- 

In the SPDC process, a pump photon (p) enters a nonlinear optical crystal 
and is annihilated, producing two daughter photons, typically called signal (s) 
and idler (i). Energy conservation dictates that uOp = To be an efficient 
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Figure 1: The dark circle illustrates the down-conversion ring. Points A, B, and 
C represent pairs of daughter photons at different azimuthal angles around the 
ring. To increase the total count rate, one should collect from multiple points 
around the ring. This idea is called spatial multiplexing. 


process, the photons must also abide by momentum conservation, a condition 
referred to as phase matching. In noncollinear SPDC, the daughter photons are 
emitted into angles on either side of the pump direction in order to satisfy this 
phase matching relation. In degenerate SPDC where ujg = this emission 
occurs in opposite pairs of directions around the pump beam forming a ring, 
called the down-conversion ring, in a plane transverse to the pump. 

In many experiments aimed at creating entangled photons pairs, the polar¬ 
ization degree-of-freedom is used because the measurement process is straight¬ 
forward and very high purity can be achieved. Because current QKD schemes 
aim at maximizing photon rates, methods of creating high brightness polarization- 
entangled sources have been investigated [sniiiaEiiiniiiii. A potential solution 
involves using Type-I phase matching, where the two daughter photons have 
the same polarization that is opposite to the pump. Recently, Rangarajan et 
al. m showed that high-brightness polarization entanglement can be achieved 
in Type-I phase matching when two thin crystals are placed together with their 
optic axes rotated 90° with respect to each other [iniiii]. Originally, the pro¬ 
cess was demonstrated in beta barium borate (BBO, a uniaxial crystal) piTO]. 
but later shown for bismuth triborate, (BiBO, a biaxial crystal) [TT] which is 
quite promising for high-brighness applications because of its higher nonlinear 
coefficient [laiis]. 

Each crystal produces the same polarization but there is essentially no which- 
path information regarding the crystal in which the photons are generated. This 
process can give rise to polarization entanglement around the entire ring, if the 
rings overlap everywhere. This enables higher achievable fluxes by allowing for 
collection of biphotons at multiple pairs of points around the down-conversion 
ring, known as spatial multiplexing and illustrated in Fig. 

To achieve entanglement around the whole ring, the collected light must be 
indistinguishable. This only occurs when the emmission patterns are completely 
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circular because any eccentricity of the ring from the first crystal will have its 
major axis perpendicular to the major axis from the ring produced in the second 
crystal. This leads to reduced spatial overlap and reduced entanglement. 

Although BiBO is a promising crystal for high-brightness SPDC applica¬ 
tions, we find that the emission pattern produced from a Type-I interaction 
with BiBO is elliptical. Eccentricity introduces a possible way of distinguish¬ 
ing photons that were born in the first crystal from those born in the second, 
thereby potentially reducing the entanglement quality around the entire ring. 

In the next section, we develop a formalism for predicting theoretically the 
emission pattern for noncollinear degenerate Type-I SPDC in BiBO and discuss 
the physical reason for elliptical emission patterns. In Sec. 3, we discuss our 
experimental setup and present our data. In Sec. 4, we discuss a theoretical 
model for predicting the eccentricity in a biaxial crystal. We show that this 
model agrees with our experimental data and that there is a wavelength that 
minimizes the eccentricity for a given set of emission angles. In Sec. 5, we 
discuss the repercussions of elliptical rings on the entanglement quality and 
brightness and further show how the spectrum of the single photons changes as 
a result of the asymmetry. 

2 Phase Matching for Type-I SPDC in BBO and 
BiBO 

The phase matching process for SPDC determines the emission direction of the 
daughter photons. In general, perfect phase matching occurs when 

kp = kg ki^ ( 1 ) 

where 

kj = ( 2 ) 

for j = where the angle- and frequency-dependent refractive index is 

given by nj{ujj^Sj), Sj is the propagation direction unit vector, and c is the 
speed of light in vacuum. We follow the geometry, notations, and conventions 
in Refs. [H ng. The crystal geometry and photon wavevectors and angles 
are shown in Fig. ga). A pump photon with wave vector kp is incident on a 
nonlinear crystal and makes angles {Op, (j)p) with the optic axis of the crystal 
denoted “OA-U” (“OA-B”) for a uniaxial (biaxial) crystal, discussed in more 
detail below. The signal (idler) photon with wavevector is emitted at local 
angles {0s{i)^ (t>s{i)) with respect to kp. For perfect phase matching, (pg = Ttt. 
The photons undergo refraction at the air-crystal interface and exit at exterior 
angles of O'^^.y 

To determine the emission angles {0s^0i,(j)s^(t)i) of the daughter photons, 
we choose frequencies {ujp,uJs,uJi) as well as the angles the pump wavevector 
makes with the optic axis {Op,(j)p). Using the Sellmeier equations for either 


3 



(a) 


OA-U 




0 60 120 180 240 300 360 

cPs (deg) 


Figure 2: (a) Crystal geometry for BBO and BiBO crystals. The pump makes 
angles {0p^(j)p) with respect to the crystal optic axis. Here, = 0. The signal 
(idler) photons are emitted at local angles (^s(i), 0s(i)) with respect to kp. Re¬ 
fraction at the interface of the crystal results in an exterior angle of outside 
the crystal, (b) The exterior emission angle versus its local azimuthal angle for 
BiBO (blue, solid line) and BBO (maroon, dashed line). For BBO, we use a 
pump cut angle of Op = 29.392° while for BiBO we use (j)p = 90°, Op = 151.563°. 
Both curves are calculated using Xp = 405 nm and = 810 nm. 


BBO or BiBO, we solve for either set of parameters by finding a solution to the 
simultaneous equation El 


{Ak^f + {Akyf + {Ak,f = 0, 


(3) 


and 

= 0, (4) 

where, for example, /Xkz = kp^ — kg^ — ki^. Evaluating Eqs. and [^involves the 
refractive index for each photon propagating through the crystal, which depends 
on the polarization of the photons, and the propagation direction. 

In a uniaxial crystal, which has a single axis of symmetry, an ordinary- 
polarized photon (o-polarized) is polarized perpendicular to the plane that con- 
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tains its propagation wavevector and the optic axis of the crystal. This photon 
experiences a refractive index that does not change with the direction of prop¬ 
agation. The extraordinary-polarized photon (e-polarized) is polarized in the 
plane of the optic axis and the propagation vector and experiences an angle- 
dependent refractive index. Type-I interactions, such as those considered here, 
include interactions where the pump is an e-polarized photon and the daughter 
photons are o-polarized, or vice versa. 

A biaxial crystal, such as BiBO, has two optic axes and reduced symmetry. 
Polarized photons are neither e-polarized or o-polarized in the sense of a uniaxial 
crystal. Depending on type of crystal and wavelength range, polarized photons 
may experience either an angle-dependent or angle-independent refractive index. 
The photons are said to be either “fast” or “slow” instead of “e” or “o” where 
fast (slow) refers to having a smaller (larger) refractive index [15]. 

Determining the fast (slow) refractive indices involves finding the length of 
the minor (major) axes of the optical indicatrix given by Fresnel’s equation of 
wave normals given by 

__^_ ^ _^__ =0 

_o/ -2 ~ _o/ -2 ' - 9 / -2 ’ y ^ 

n ^{uj^s)—nx n ^[uj,s)—ny n ^[uj^s)—nz 

where Ux^riy, Uz are the refractive indices in each principle direction of the crys¬ 
tal at a given vacuum frequency and n{uj,s) is the refractive index in a given 
direction with unit vector s. For a negative biaxial crystal, such as BiBO, 
'^x < '^y < Solving Eq. for n{uj^s) using the approach described in 
Ref. [H], we obtain two solutions, one for each polarization (fast and slow). In 
BiBO we find that the pump photons experience the fast refractive index that is 
angle-independent for wavelengths in the UV and blue part of the spectrum. In 
contrast, the daughter photons in the red and NIR part of the spectrum experi¬ 
ence the slow refractive index that is angle-dependent. In BBO, pump photons 
in the blue part of the spectrum travel as e-polarized photons, and experience 
an angle-dependent refractive index, while the signal and idler photons are o- 
polarized and experience an angle-independent refractive index. This effect is 
not due to the uniaxial versus biaxial nature of the crystals, but simply due to 
the angle dependence on the refractive index over certain wavelength ranges for 
each particular crystal. 

Our analysis of the elliptical emission pattern for BiBO has not been noted or 
observed previously. This may be due to the fact that most previous quantum 
optics experiments have been conducted with BBO, which produces circular 
rings as observed and predicted below. This is not the case in BiBO due to 
the angle-dependent refractive index experienced by the daughter photons. In 
Fig. I^b) we plot the external angle as a function of its azimuthal angle (j)s 
for BiBO (blue, solid) and BBO (maroon, dashed). For BiBO, we observe that 
6>' varies with (j)s so that 6>' has a larger value at (j)s = 0° and 180° than in the 
(j)s = 90° and 270° directions. This leads to an elliptical emission pattern for 
this crystal. For BBO, 0'^ is a constant as a function of which implies a 
circular emission pattern. In our analysis, we ignore birefringent walk-off of the 


5 





beams due to its negligible contribution for the thin crystals considered here, as 
discussed in detail in the Appendix. 


3 Experimental Results 

We image the spatial intensity patterns from both BiBO and BBO using the 
experimental setup shown in Fig. A 405-nm-wavelength, continuous-wave 
laser pumps either a BBO or BiBO crystal. Both types of crystals are thin 
(0.8 mm) in comparison to the Rayleigh length and compared to the transverse 
extent of the imaged rings - a requirement for polarization entanglement. As 
the crystals become thick, which-path information begins to degrade the entan¬ 
glement because the rings become distinguishable. We test two different BiBO 
crystals each with a different set of crystal cut angles {0p^(j)p) designed to be 
phase matched for \p = 405 nm and Ag = = 810 nm with an exterior open¬ 

ing angle of ^ 3°. These two crystal cut angles are the angles for which it is 
straight-forward to create an optical beam that propagates through the crystal 
with negligible walk-off. That is, <pp is chosen so that the pump, signal, and 
idler photons essentially propagate as a slow or a fast wave. We can further tilt 
the crystal around Op to tune the opening angle. One crystal has cut angles of 
(^p = 90°,Op = 151.7° while the other has (pp = 0°,0p = 51°. The BBO crystal 
has a pump cut angle of Op = 29.3°. 

Object plane 



Figure 3: Experimental setup. A 405 nm laser (Omicron, 

LDM405.120.CWA.L.WS, < 0.02nm bandwidth FWHM) pumps either a 
BiBO or BBO crystal (Newlight Photonics). The down-converted light exits 
the crystal and propagates 12.4 cm before it passes through a 40 cm focal 
length lens (Thorlabs LAC726B). This lens functions to direct the rays from 
the down-conversion ring to an object plane that we image onto the camera. 
We image the ring ^ 10 cm after this lens by placing a zoom lens (Navitar 
Zoom 7000E) on an EMCCD (Andor iXon^^) ^ 1 m away from the lens. The 
magnification of this imaging system is 8.6. Each pixel on the camera chip 
is 24 /im X 24 /im, the sensitive area of the chip is 3 mm x 3 mm, and we 
cool the chip down to —70°C to reduce dark noise. We put a 10 nm bandpass 
optical filter (Andover 810ES10-50) before the camera lens to select only nearly 
degenerate wavelengths. 
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Making precise measurements of the eccentricity requires having the well- 
calibrated imaging system, depicted in Fig. To achieve this, we perform 
the following procedure: first we remove the crystal and the 40 cm focal length 
lens, leaving only the laser, a set of apertures for alignment, steering mirrors, 
and camera-lens system. We then place a flat mirror against the camera lens to 
ensure the back-reflected light is going straight back through the apertures. We 
then place a target in the object plane. The target is a flat piece of metal with 
a 20-mm-diameter ring scored in the surface of the metal, where the diameter 
tolerances is ~ 25 fim. This machined ring has a hole in the center so that 
we can easily align it with the laser beam path and check its back reflections 
for tilt. We check the eccentricity of the target and make small adjustments to 
the camera’s position and tilt until we minimize eccentricity. The eccentricity is 
caused primarily by any amount of tilt in the system, which arises from imperfect 
alignment. Astigmatism and coma also arise from an imperfectly aligned optical 
system, although these are negligible compared to the tilt for a well-aligned 
optical system. We then replace the crystal and check back reflections with a 
mirror on the camera lens. Finally, we add in the lens and ensure its alignment 
by checking back reflections. 

We collect images for each crystal and fit the observed emission pattern 
with an elliptical function in two transverse dimensions with a Gaussian profile 
in the longitudinal dimension. The free parameters of our model include the 
height of the Gaussian peak, background counts (offset of the Gaussian from 
0), major/minor axes for ellipse, width of Gaussian and location of the center 
point. We calculate the eccentricity of the ellipse by 



where b (a) is the major (minor) axis of the ellipse. 

Figure shows the transverse spatial intensity pattern for down-conversion 
rings from both BBO and BiBO crystals. For BBO, the intensity pattern of 
the down-conversion ring is essentially circular (Fig. [^(a)), e = 0T0.013 where 
the error is a combined statistical and systematic error of 0.013, which will be 
discussed below. Hence, our results indicate the pattern for BBO is circular 
to within our experimental uncertainties. The range of opening angles that 
are phase matched for BBO is smaller than that for BiBO because the slope 
of the opening angle versus wavelength is larger in BBO due to the slope of 
the refractive index versus angle being steeper. This leads to a thinner down- 
conversion ring for BBO. For BiBO cut at phase matching angles (pp = 90°, Op = 
151.7°, (Fig. (b)), the eccentricity is greater than that for BBO, with e = 
0.172 ± 0.019. For BiBO, with Pp = 0°, Op = 51° (Fig. [d] (c)), e = 0.367 ± 0.012 
and is easily seen by eye. 

Table 1 gives the results of our observations and analysis. We calculate the 
major/minor axes using the experimental values for 0'g{(j)s = 0)/6>'(0^ = 90) and 
Eq. Using the procedure outlined in Sec. 2, we calculate theoretical values 
for the eccentricity using an exterior angle for either the major or minor axis 
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Figure 4: SPDC emission pattern for BBO and BiBO. The emission pattern 
from the camera is plotted versus the two transverse dimensions {x and y) 
where y is in the direction perpendicular to the optical table. The patterns 
from the camera have been scaled up by the magnification and converted from 
pixels into cm. The red, solid lines are ellipses with the major and minor axes 
taken from the fit parameters, (a) Down-conversion ring from BBO crystal is 
nearly circular with an eccentricity of 0.013. (b) Down-conversion ring BiBO 
crystal cut at phase matching angle of (Op = 151.7°, = 90°) has a higher 

eccentricity with the major axis in the x-direction. (c) Down-conversion ring in 
BiBO crystal cut at phase matching angle of (Op = 51°, = 0°) has a large 

eccentricity. 


given in the first column. The exterior angle is what we measure experimentally 
and is simply the opening angles in the crystal propagated outside the crystal 
Snell’s law. We use this set exterior angle to calculate the Op and (pp that gives 
this value at either = 0 or = 90°. In this case, instead of knowing Op and 
calculating Os and we know the latter and calculate the former. Measuring 
the external emission angle is more straightforward than measuring the angle 
the pump beam makes with the crystal optic axis. Once we determine Op^ we 
use it to predict the emission angle around the entire ring. Using Snell’s law, 
we propagate these angles outside the crystal and find the eccentricity through 
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Table 1: Experimental Data 



BBO 

BiBO 
(/)p = 90° 

BiBO 

</>p = 0° 

Experiment (°) 

e 

4.12358 

4.12382 

0.013 

4.0294 

4.10893 

0.172 

4.05449 

4.31798 

0.360 

Theory e 

0 

0.166 

0.361 

Error 




Statistical 

0.011 

0.001 

0.001 

Systematic 

0.007 

0.019 

0.012 


the relation, 


/ = 0 °)] \ 
Vtan[6»'(?!'s = 90°)]/ 


(7) 


where the ± accounts for the possibility that the axes at = 0 may be the 
major or the minor axis. We use Eq. to calculate both experimental and 
theoretical values for e. Eor the experimental values, we have measured both 
tan[6>'(05 = 90°)] and tan[6>'((/)s = 0°)], while for the theoretical values, we use 
only one of these, either tan[^'(0^ = 90°)] or tan[^'(0^ = 0°)] and calculate the 
other by the procedure described above. The error in our experiments includes 
a statistical error from the fitting process and a systematic error, which is most 
likely due to small astigmatism in the system. We calculate the systematic error 
by taking repeated measurements at different times which requires realignment 
using the alignment procedure outlined above. 


4 Theoretical Model 

In this section we outline a theoretical model for predicting the approximate 
eccentricity of the emission pattern in each crystal. Because our emission an¬ 
gles are small {Og^Oi « 1), a small-angle approximation around the collinear 
case {Os = = 0) is a natural method for obtaining an approximate analytic 

solution. This allows us to gain better understanding of the parameters that 
affect the eccentricity. Starting from the phase matching equations in the y and 
2 : directions, we find that the phase matching occurs when. 


cos(6>s) + cos(6>i) — 2np = 0, (8) 

Us sm{0s) + rii sm{0i) = 0. (9) 

We expand sin(6>j), cos(6>j), n^, and rii to second order in Oj. The approxi- 
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mate expressions of the refractive indices are given by, 



( 10 ) 



( 11 ) 


where Op = Op^ + SOp is the pump phase matching angle, Op^ is the pump 
phase matching angle for the collinear degenerate case and SOp is a small vari¬ 
ation around that angle. The notation Ug describes the refractive index for the 
collinear degenerate case, and Ug has replaced everywhere because for the 
collinear degenerate case, they are equal. The change in refractive index with 
angle is calculated from the expression for the frequency and directional depen¬ 
dent refractive index. These quantities indicate how quickly the refractive index 
changes with the various angles. Inserting Eqs. pTjand ?? and the trigonometric 
functions into Eqs. |^and[^ we solve for angles Og and Oi. We are interested in 
solving for SOg and SOi which are small changes in the emissions angles around 
the values of 0g{(l)g = 90°) and 0g{(l)g = 0°), which are the angles corresponding 
to the directions of the major and minor axes. We use the expression 




( 12 ) 


to solve for S0s{(j)g = 0°) and S0g{(j)g = 180°). We calculate the eccentricity 
inside the crystal to be 



(13) 


The eccentricity arrises from the photons experiencing the angle-dependent re¬ 
fractive index, so the eccentricity in the emission pattern is only due to the 
eccentricity inside the crystal. Once the photons propagate into free-space, 
they experience no angle-dependence and simply follow the trajectories of the 
external angles dictated by Snell’s law. Therefore, the eccentricity inside the 
crystal is the only contribution to the eccentricity experimentally measured. 

Using Eq. we determine the predicted value of the eccentricity using the 
measured values for SOp. These are shown in the final columns of Table 1. Eor 
small values of SOp^ we also express the eccentricity in terms of the derivatives 
of refractive indices as. 



( 14 ) 
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Figure 5: Plot of each term in Eqj^ versus wavelength, “1” (blue, short-dash) 
{d‘^ns/dOg)^^^^ “2” (red, long-dash) ^ “3” (green, solid) 2fns^ 

[{dns/dOs)^^^^]‘^. For this plot we use (pp = 90°, Op = 152.077°. 


where the superscripts of denote the azimuthal angle at which the terms are 
evaluated. We plot the three terms in this approximation in Fig. For this 
plot, the difference between the curves “1” and “2” is the dominant contribution 
to the eccentricity. Where these curves intersect is where the eccentricity is at a 
minimum in Fig. This expression allows us to determine why eccentricity is 
larger for certain crystal cuts than others. We examine the relative size of each 
of the three terms in Eq. 14 for the case where = 0 and (j)p = 90°. We find 


that the difference of the first two terms in Eq. 14 for (pp = 0 is very large in 
comparison to the difference of the first two terms for the pp = 90° case. This 
ultimately leads to the larger eccentricity we observe for the (pp = 0 crystal cut 
in BiBO. 

Using this method, we also calculate the value of the degenerate wavelength 
(As = \i) that minimizes the eccentricity inside the crystal. Again, minimizing 
the internal eccentricity will minimize the eccentricity measured in the external 
emission pattern. Each of the parameters in Eq. is wavelength-dependent 
because the refractive indices for the daughter photons in BiBO are angle- 
dependent. We find a minimum wavelength by performing a root-finding routine 
on Eq. [^with wavelength as the independent parameter. We also find the ec¬ 
centricity as a function of wavelength plotted in Eig. [^for three different values 
of Op. This plot shows that there is a wavelength for each pump tilt angle for 
which the eccentricity can be minimized. According to the plot, the eccentricity 
becomes large away from the 710-780 nm range. Eor our degenerate wavelength 
of 810 nm, this figure agrees with the value of the eccentricity we observe. This 
plot indicates which wavelengths are best at minimizing eccentricity in a given 
crystal cut for Type-I degenerate SPEC. 
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Figure 6: Eccentricity versus wavelength for three different values of Op^ “1” 
(blue, short-dash) Op = 152.071° “2” (green,solid) Op = 151.378°, “3” (red, 
long-dash) Op = 149.21°. For all three plots there is a wavelength for which the 
eccentricity of the down-conversion ring is minimized. We show our experimen¬ 
tal data for a pump wavelength of 405 nm and a degenerate down-converted 
wavelength of 810 nm (purple) with associated error bars. 

5 Impact on entanglement and count rate 

One potential drawback of an elliptical emission pattern is that it may decrease 
the entanglement quality and/or the entangled photon count rate. Because the 
rings are elliptical, they no longer perfectly spatially overlap around the entire 
ring, leading to potential degradation in entanglement-purity, and limiting the 
ability to multiplex many channels around the down-conversion ring. Lack of 
spatial overlap is the most obvious way in which which-path information is 
revealed. The elliptical shape may also cause a change in the joint spectrum 
for the photons, which leads to a degradation in entanglement purity. We note 
that spatial overlap can be corrected for by collecting light in between the ring 
centroids at the cost of count rate. However, at this location, the spectra from 
the two rings is different, leading to decreased entanglement purity. We study 
each of these effects in different ways where we assume that the biphotons are 
coupled into single mode fibers. 

To estimate whether the entanglement quality is degraded, we determine 
the spectrum of the down-converted light for each crystal for single mode fiber 
collection. Using single mode fibers is important because collection of a single 
spatial mode for either the signal or idler photon means that the twin photon 
should also be projected onto a single spatial mode. The entanglement purity 
can then be determined from the overlap integral of the spectral intensity of 
these two single modes. Calculating the spectra from the midpoint between the 
rings at two locations and computing the overlap integral provides a measure 
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of distinguishability. A smaller overlap integral represents lower entanglement 
quality because the difference in spectra between the crystals means it is possible 
to distinguish which crystal the photons came from. 

Using the formalism developed in our previous work HZ], we determine both 
the joint spectrum for the biphoton wavefunction as well as the singles spectrum 
for each signal and idler, where we assume single mode fiber collection. Figure 
(b) shows the joint spectrum at location B in Fig. (a). Here the overlap is 
minimal for each crystal in the pair when the crystals are cut at (pp = 90°, Op = 
151.7° (solid lines) and pp = 0,6>p = 51° (dashed lines). In these calculations 
we choose one of these points, point A (Fig. (a)) to have maximal overlap 
between the two rings, leaving point B to have minimal overlap. We choose 
the minimal-overlapped case because the entanglement purity is lowest due to 
the mismatched overlap. For the minimal overlap case, we choose the collection 
mode to be in between the two rings so that we collect equal count rates from 
each and so that the spectra from the two will be most similar. We find that, 
for both BiBO crystal cuts, the spectra are similar for the two crystals, and the 
entanglement quality is not substantially degraded. The spectra at point A are 
identical, while at point B we find the overlap integral is 99.4-99.97% for each of 
the two crystal cuts. The overlap integral for two entangled states is a measure 
of the entanglement purity of the states involved. 

Even though the eccentricity of the ring does not affect the entanglement 
purity, it does affect the count rates. Integrating the spectra in Fig. (b) over 
Acj gives the total joint count rate. For the crystal cut at pp = 0, = 51°, the 

spectra has a much lower amplitude and integrating over any band of frequencies 
gives a lower count rate compared to the crystal cut at pp = 90°,6>p = 151.7°. 
Higher eccentricity causes the rings to be more spatially separated in the former 
case, resulting in the lower joint count rate. Additionally, the crystal cut for 
this high eccentricity {pp = 0,6^^ = 51°) has a lower effective nonlinearity (deff) 
resulting in an even lower count rate so that this crystal cut is clearly not opti¬ 
mal for high-brightness applications. This may not be the case for other crystals 
however. For SPDC applications in BiBO that require both a high count rate 
as well as high purity entanglement, we find that the eccentricity does not de¬ 
crease the entanglement purity significantly, however, the count rate is lower for 
collection of the midpoints between emissions patterns with greater eccentricity. 
Applications aiming for ultra-high entanglement purity (> 99.99%) are limited 
to the locations around the ring that can be made to overlap completely. 

6 Conclusions 

In conclusion, we show both experimentally and theoretically that the emis¬ 
sion patterns from down-conversion in BiBO have an elliptical shape, while 
the pattern from BBO has a circular profile. We show that this difference in 
shape of the down-conversion rings depends on whether or not the daughter 
photons experience an angle-independent or angle-dependent refractive index. 
Although, in hindsight, the results we present may seems obvious, this is to our 
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Figure 7: (a) Points A and B represent the points of maximal and minimal 
overlap respectively. Red inner ring lines and blue outer ring lines illustrate the 
center of the emission patterns from each crystal in the crystal pair, where we 
have ignored the finite width of the actual ring. The eccentricity of these rings 
has been exaggerated for illustration purposes. Left figure depicts crystal cut at 
(f)p = 0,0p = 51° and right figure depicts crystal cut at (pp = 90°, Op = 151.7° (b) 
Joint spectra for each crystal cut at location B for each crystal cut. The y axis 
on the plot is the differential rate per frequency and the x axis is an angular 
frequency shift around the degenerate frequency. Solid lines (green (2) and gold 
(1)) represent the joint spectral for each crystal when the crystals are cut at 
(j)p = 90°, = 151.7° while the dashed lines (blue (3) and purple (4)) represent 

the joint spectra from each crystal when the crystals are cut at pp = 0, Op = 
51°. The pump wavelength for this simulation is 405 nm and the degenerate 
wavelengths 810 nm. The vertical dashed green lines in the plot indicate a 20 
nm bandwidth for reference. 


knowledge, the first time this effect has ever been characterized and its poten¬ 
tial impact discussed. We also present a theoretical method for calculating the 
eccentricity of the down-conversion ring for biaxial crystals find that there is an 
optimal wavelength for which the eccentricity is a minimum and close to zero. 
We demonstrate that the elliptical nature of the rings does not reduce entan¬ 
glement purity, but reduces joint counts rates significantly for down-conversion 
patterns with a larger eccentricity. 
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Appendix 

Birefringent walk-off occurs in a medium when, for a particular polarization, the 
momentum vector and the Poynting vector separate from each other. Walk-off 
occurs in both uniaxial and biaxial crystals, although calculating the walk-off 
angle in the biaxial case is more challenging due to the reduced symmetry of 
the crystal. We calculate the walk-off angles using the method in Ref. m, and 
find the contribution to the eccentricity from walk-off is negligible, due to the 
fact that the crystals considered here are so thin. The crystal thinness means 
that photons that are walking off different amounts in different directions will 
not cause a large change in eccentricity because their propagation distance to 
the crystal face is quite small. 

In BiBO, the daughter photons experience an angle-dependent refractive 
index, and will therefore have Poynting vector walk-off. Using Eqs. (20)-(22) in 
Ref. [15] we determine the Poynting unit vector N and unit propagation vector 
k. We use these to determine the walk-off angle as a function of azimuthal 
angle around the original pump direction. We find walk-off angles for the 
4>p = 90°, Op = 151.56° crystal cut at = 0 and = 180° are 3.19° and 
3.51°, respectively. The walk-off angles for (jyg = 90° and = 270° are both 
3.36°. We map the Poynting vectors and momentum vectors onto the crystal 
exit face and compare their eccentricity. We find that the Poynting vectors map 
a ring with eccentricity of 0.1680, while the momentum vectors map a ring with 
eccentricity 0.1685, which is a 0.3% difference. 

In free space, the Poynting vector and momentum vector must be parallel, 
so the trajectories of these vectors to the image plane are identical. As we have 
shown in the previous paragraph, the eccentricity of the ring as it is about to 
exit the crystal is negligible. The photon path outside the crystal is much larger 
than the one inside the crystal and as the Poynting vector and momentum vector 
are parallel in free-space, the dominant contribution to the eccentricity comes 
from the paths taken outside the crystal. Essentially, the overall impact of the 
walk-off from intracrystal angles is very small due to the fact that the crystal 
length is much smaller than the crystal - lens distance. 
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